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ABSTRACT 

A thermodynamic foundation using the concept of internal stale variables is given for a 
general theory of viscoplasticity, as it applies to initially isotropic materials. Three funda- 
mental internal state variables are admitted; they are: a tensor- valued back stress for 

kinematic effects, and the scalar-valued drag and yield strengths for isotropic effects. All 
three are considered to phenomenologically evolve according to competitive processes 
between strain hardening, strain-induced dynamic recovery, and time-induced static recovery. 
Within this phenomenological framework, a thermodynamically admissible set of evolution 
equations is put forth. This theory allows each of the three fundamental internal variables to 
be composed as a sum of independently evolving constituents. 


1. INTRODUCTION 

The theoretical development of viscoplasticity has its origin with the works of BINGHAM & green [1], 
HOHENEMSER & PRAGER [2], OLDROYD [3], ODQVIST [4], STOWELL [5], and PRAGER [6], whose models do not 
contain evolving internal state variables. The field started to gain inertia in the mid 1960’s when internal state 
variable models began to appear in the theories of perzyna [7] and Armstrong & Frederick [8]. With the 
advent of the computer, rapid advances were made in the 1970’s with the modeling efforts of BODNER & partom 
[9], HART [10], MILLER [11], PONTER & LECKIE [12], CHABOCHE [13], KRIEG et al [14], and ROBINSON [15]. 
Theoretical refinements have continued to occur throughout the 1980’s in the models of chaboche et al. [16], 
STOUFFER & BODNER [17], MARQUIS [18], VALANIS [19], WALKER [20], SCHMIDT & MILLER [21], CHABOCHE & 
ROUSSEUER [22], EISENBERG & YEN [23], BROWNING et d . [24], ESTRIN & MECKING [25], KREMPL et al. [26], 
LOWE & MILLER [27], ANAND & BROWN [28], HENSHALL et al. [29], MCDOWELL & MOOSBRUGGER [30], NOUAILHAS 
[31], and freed & walker [32]. Reviews on various aspects of viscoplasticity have been written by perzyna 

t Prepared for submission to the ASME Journal of Applied Mechanics. 
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[33], WALKER [20], CHAN et al [34], LEMAITRE & CHABOCHE [35], SWEARENGEN & HOLBROOK [36], EFI7S & 
JOLLES [37], miller [38], and chaboche [39]. Although this listing is by no means complete, it does provide 
the reader with a represen ti live bibliography of the work done in the field of viscoplasticity for initially isotropic 
materials. 

This paper presents the development of a general, nonisothermal, viscoplastic theory from physical and 
thermodynamical considerations. The basic constitutive equations, evolution equations, and dissipation inequali- 
ties are derived for a class of materials whose virgin state is stress free and isotropic. Specific functional forms 
for the evolution equations are then proposed as they apply to a theory of viscoplasticity. Special attention is 
given to the intrinsic dissipation inequality. Here three fundamental internal state variables are considered. One 
accounts for kinematic or flow-induced anisotropic effects (i.e., the tensor-valued back stress), while the other 
two account for isotropic effects (i.e., the scalar-valued drag and yield strengths). Each is considered to evolve 
according to competitive phenomenological processes between strain hardening, strain-induced dynamic 
recovery, and time-induced thermal recovery, in accordance wilh current philosophies in the field ( cf review 
articles listed above). Of particular interest is the fact that each of the three basic internal variables can be com- 
posed as a sum of separately evolving constitutents, which are the primary internal variables. Aspects of this 
feature have been used by chaboche et a /. [16,22,35,39], miller et al [21,27,29,38], and others in their visco- 
plastic models. 

We shall consider the thermo-mechanical behavior of viscoplastic materials; in particular, polycrystallinc 
metals and their alloys. Each element of the material is assumed to be isotropic and to carry no stress in its ini- 
tial virgin stale. As the material deforms, however, anisotropies may become induced. Only kinematic anisotropy 
is considered in this paper. Small material displacements and rotations shall be considered to make up the defor- 
mation of the material. In a Cartesian reference configuration, the strain e tJ is taken to be composed of elastic 
zfj (reversible - includes thermal strain) and inelastic or plastic zfj (irreversible) parts such that 

By = By + e* 0) 

and there is no inelastic strain in the stress-free virgin state. 

2. THERMODYNAMICS 

The changing internal structures of a material element are characterized by its stale {z‘j,a^T t V t T}. The 
clastic strain is a measure for elastic changes in the internal structure. The internal state variables 
a ^ (^=1,2, . . . ,n) are measures for inelastic changes in the internal structure. The temperature T is a meas- 
ure for the heat in a material element. And the thermal gradient V,T is a measure for heat fluxing from the 
clement. An inelastic material response brought about by a change in the internal variables a ^ affects the 
future response of the material. If the a ^ are taken to be scalar valued, then nonintersecting subsets of the a $ 
may be taken to be components of higher-order, irreducible, even-rank tensors associated with a reference 
configuration [40]. 
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For any quasistatic adiabatic process (which is also a thermostatic or reversible process), Gibbs equation is 
given by 

dS = j(du- Gjjde.j) (2) 

where S is the entropy, u is the internal energy, c is the applied stress, and the set } defines the 

slate. Entropy and temperature both exist in thermostatics, because dS is an exact differential and T > 0 is 
its associated integrating factor [41], It is from this relationship that one derives the theory of thcrmoelasticity. 
Here repeated Latin indicies are summed from 1 to 3 in the usual manner. 

For any thermodynamic process, be it reversible or irreversible, the Gibbs equation of thermostalics gen- 
eralizes to the inequality 

S > y(ii - Cijkij + (3) 

where q, is the heat flux, and the set {5',e t y,V J T} defines the state. This relationship follows from the conser- 
vation of energy and the Clausius-Duhem inequality [42]. Here the existence of both entropy and temperature 
must be assumed to continue into the domain of irreversible processes, as rigorous justification for their 
existence continues to elude theoretical proof [43]. It is from this relationship that we derive our theory of 
viscoplasticity. 

Thermostatics (equation (2)) addresses reversible processes where material response depends only on the 
initial and final states. Thermodynamics (equation (3)), on the other hand, addresses irreversible processes 
where material response not only depends on the initial and final states, but it also depends on the path traversed 
in state space. 

Introducing the thermodynamic potential known as the Helmholtz free energy 

V - u - ST (4) 

into equation (3), and using the decomposition of strain given in equation (1), one obtains the inequality 

Oij ifj > y + ST - CTy i-j + jq, V, T (5) 


From our experience with thermostatics, this relationship suggests that the set (e£,e£,7\V- T) should be the set 
of independent state variables. However, experimental evidence dictates that e* is not an admissible indepen- 
dent state variable (cf. onat & fardshishhh [44]), Rather, one ought to introduce a set of internal state vari- 
ables a ^ (£=1,2 n) such that the set {e^,a^,T,V-r} is taken as the set of independent state variables 

[45]. Therefore, considering the potential 


y = \j/(e<,a f : 1 7',V l T) 


( 6 ) 


equation (5) becomes 


Gij&fj — 


’iSL - 
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But because this inequality must be valid for arbitrary changes in the observable variables £,*, T , and dV^T/dt , 
as the contributions from their associated terms are reversible, one obtains the following system of equations 
governing any thermodynamic process. They are [45]: the constitutive equations 




dr ’ dv.r 


( 8 ) 


the dissipation inequality 1 * 

~ y?. V . T * 0 

5=1 1 


(9) 


and the evolution equations 

= efj(e& 0 v Ty k T) , = d^tuAyJ k T) , q t = ? I (e£,a T T,V t 7') (10) 

where /4^ (5=1,2 n) are the thermodynamic forces or affinities conjugate to the thermodynamic displace- 

ments a $ or fluxes a^. Notice that the Helmholtz free energy is independent of the thermal gradient [42,45]. 

Our past observations of material response allow us to consider q t as being independent of and 

T, explicitly, and to also consider both if’ and as being independent of V, T, explicitly. Hence, the total 
dissipation inequality given in equation (9) can be separated into the intrinsic dissipation inequality given by 

Vijifj-YAiP K~° (11) 


and the thermal dissipation inequality given by 

yft V.T <0 


( 12 ) 


U suall y, the evolution of heat flux is represented by Fourier’s equation 

<7, = -kijVjT (13) 

where the symmetric tensor for thermal conductivity must be positive definite in order to satisfy the thermal 
dissipation inequality. 


3, CONSTITUTIVE EQUATIONS 

The constitutive equations derived from thermodynamic considerations in the previous section are very 
general. However, by choosing a particular form for the thermodynamic potential y based upon oiir past 
experiences, these general relationships become specific and useful. 


1 Some authors (eg., COLEMAN A GURTIN [42]) would divide this inequality by temperature and set it equal to a 

non-negative term they call the rate of entropy production, but this is an unnecessary definition in our development. 
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Most existing viscoplastic models (without damage) belong to the class of materials characterized by the 
thermodynamic potential 

V - = V&jJ) (14) 

such that 

V* = ki ~ ZijDijki a kikT - C £ T In 

/ = -° (15) 

where T 0 is the initial value of temperature with AT = T-T 0 , C E = T dS (z t j t T)ldT is the heat capacity at 
constant elastic strain, a tJ are the thermal expansion coefficients, D ijkl are the elastic moduli, and is the 
symmetric matrix of hardening coefficients (which is taken to be positive definite). These material constants 
typically vary with temperature to some degree. It is in this sense that the plastic potential depends on 
temperature. This is an implicit dependence rather than an explicit one, however, as we imply by the notation 
(;7). The fact that the plastic potential is taken as a positive-definite quadradic form is a stronger require- 
ment than is actually necessary, since only the property of convexity is required of this potential. 

By substituting the thermodynamic potentials of equation (14) into the constitutive equations (8), and 
using the decomposition of strain given in equation (1), one obtains the following system of equations which 
govern any process for the class of viscoplastic materials considered. They are: the constitutive equations 

£47 = DijkiGu + o.ijAT + efj 



S = a ijCij + C a In 



( 16 ) 


2 Hfyr a y 

the dissipation inequalities 

V * T ^ 
5=i 


and the evolution equations 

O’ = ifauAyJ) , d (t = d (> (a k ,A v T) , d. = < 18 ) 

where C a = T 35 (c i; ,T)/dT, and therefore C° = C E + a t) D l]kl a u in the limit as AT 0, which denotes the 
heal capacity at constant stress or pressure. Here it is assumed that 5=0 when AT = a,,- = 0. These govern- 
ing equations reduce to those of thermoelasticity [46] in the absence of inelastic deformation. The thermo- 
dynamic forces replace the thermodynamic displacements as the arguments of the evolution equations without 
loss in generality because the constitutive equations for internal state are linear. This fundamental set of 
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equations forms the basis upon which we develop our viscoplastic theory in the next section. 

Since this paper is directed towards materials that are isotropic ip their initial virgin state, one requires 
that 

°i]ki = “8,* 5,, - jb ^ , a ij = abij , k u = kb u 

A ^ = Hty a y — » A = ha for all scalar state variables a y -> a (19) 

t+6 

f hy a y = H a u = Ha x) for all tensor state variables a y -» a i; = a Jt , a u = 0 

Furthermore, isotropy also requires that the inelastic strain, and all internal state variables that are associated 
with second-rank tensors, to be zero-valued in the initial virgin state, i.e., 

e5(/=Q) = 0 , a i} {t = 0) = 0 (20) 

Although internal state variables associated with fourth and higher even-rank tensors are admissible, they arc not 
required in this work. Here E is the elastic modulus, v is the Poisson ratio, a is the coefficient of thermal 
expansion, and k is the thermal conductivity, all of which vary with temperature, typically. The hardening 
coefficients for the internal variables h and // may vary with temperature, too. The Kronecker delta 5 V is 
defined as 1 whenever i=j , otherwise it is 0 valued. Since all internal state variables that are second-rank 
tensors must be irreducible [40], they must either be symmetric and traceless (i.e. 9 deviatoric) or antisymmetric, 
but only deviatoric tensors shall be considered here. 

4, VISCOPLASTIC THEORY 

In the preceeding section, structure was given to the constitutive equations. In this section, structure is 
given to the evolution equations resulting in a general theory of viscoplasticity. Once again wc use our past 
experiences to guide us, this time for the purpose of proposing functional forms for the evolution equations. 
One docs not have complete freedom, however, in choosing these forms. They must, by necessity, satisfy the 
thermodynamic constraint of intrinsic dissipation (see equation (11)); just as Fourier’s equation (13) must satisfy 
the thermodynamic constraint of thermal dissipation (which constrains the thermal conductivity k i; to bo posi- 
tive definite). 

The total dissipation inequality of equation (9) has the form 

^>0 (21) 

5=1 

where the X 5 represent the thermodynamic affinities {a^-Ay-V^/T} and the represent their associated 
thermodynamic fluxes characterized through the evolution equations (18). As the simplest case, the 

linear evolution equations of onsager [47] 

y (22) 

Y=1 
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apply to some neighborhood - be it large or small - of every equilibrium state, and as such, the total dissipation 
of the material element is a minimum in that neighborhood [48]. Here the symmetric matrix of coefficients L ^ 
must be positive definite in order to satisfy the thermodynamic dissipation inequality given in equation (21). 
Fourier’s evolution equation (13) for heat flux satisfies such a linear relationship; however, a linear relationship 
for the evolution of inelastic strain only exits for Coble and Nabarro-Herring creep in metallic materials (cf. 
ASHBY [49]). Neither of these creep mechanisms are of interest to viscoplasticians, because the stress levels 
associated with Coble and Nabarro-Herring creep are too small for the viscoplasticians’ intended applications. 
For those applications of interest to the viscoplastic modeler, the inelastic strain rate is a highly nonlinear func- 
tion of stress. Hence, the thermodynamics of Onsager does not apply to viscoplasticity. We are therefore left to 
consider a more phenomenological approach. 

BRIDGMAN’S [50] experimental results demonstrate that the evolution of inelastic strain in metals is, to a 
good approximation, insensitive to moderate levels of hydrostatic pressure. Consequently, the evolution equa- 
tions (18) for inelastic strain and the internal state variables take on the more restricted form 

Hi = tf(S k iAy,T) , a* = d^(S kl A v T) (23) 

where S tJ = < 5 tj - (a u /3)5 i y is the dcviatoric stress (not to be confused with the entropy S). 

The stress dependence of inelastic flow is strongly influenced by the material’s anisotropy, which may be 
either inherent, or flow-induced, or both. In this paper only flow-induced anisotropy is considered, and it is 
introduced after the manner of PRAGER [51]. This is done through an internal variable B i} called the back or 
internal stress, which is a symmetric and traceless tensor ( i.e ., deviatoric), and therefore irreducible in accor- 
dance with the theoretical results from group theory of GEARY & ONAT [40], To assure an initial condition that 
is isotropic, 0) = 0. 

von MISES [52] introduced the concept of a plastic potential F in the flow equation for inelastic strain, 
which we write as 


<24> 

where dFldOij defines the direction (in unit length) of inelastic straining, with \e?\ providing its magnitude. 
The existence of this potential followed naturally from von MISES’ [53] prior definition of yield, which is 
equivalent to introducing a second-invariant (or quadradic) norm that we write as 

n = V2/3 ifjtfj , M\ = V3/2 J.jJ.j (25) 


where the tensor J^e [S i; £ tJ -B^ ) . More recently, rice [54] demonstrated the physical existence of such a 
potential function based upon the mechanics and thermodynamics of dislocation slip. Most viscoplastic models 
take F =|S -B | such that equation (24) becomes 


■,r. = B,J 

2 \S-B\ 


(26) 
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which is compatible with the kinematic constructs proposed by PRAGER [51] in his plasticity model. The choice 
of this plastic potential provides a reasonable approximation to the actual shapes of observed yield and flow sur- 
faces [55]. The norms given in equation (25) are scaled for tension. They could have just as easily been scaled 
for shear by replacing the coefficients that appear under the radical signs with a 1/2 for S tJ £ tJ , and S t) B,, , 
and a 2 for ejy. 

ZENER & HOLLOMON's [56] experimental results demonstrate that the functional dependence for the magni- 
tude of inelastic strain rate |e p | can, to a good approximation, be decomposed into the product of functions 

|e p (Sij A P ,T) I = 0(n Z (Sij A P ;T) > 0 (27) 

where 0 > 0 is the thermal diffusivity, and Z > 0 shall be called the Zener parameter. The thermal diffusivity 
is often represented as an Arrhenius function of temperature, which is valid over a significant but specified tem- 
perature range (cf. miller [11]). The functional dependence of the Zener parameter is of particular interest to us 
in this paper. This parameter is a temperature normalized function for the magnitude of inelastic strain rate that, 
at most, has only an implicit temperature dependence. 

In addition to the kinematic variable, or back stress B tJ , there arc two isotropic variables introduced into 
the general structure of our theory; they are the yield strength Y > 0 and the drag strength D > 0 (not to be 
confused with the elastic moduli D ijkt ). These three internal variables interact with the deviatoric stress 5,, in 
such a way that the Zener parameter of equation (27) is considered to have the following functional dependence 


in accordance with existing viscoplastic models. Here a 1 ' = \S-B\-Y is the viscous stress governing the ine- 
lastic material response, with Z(0) = 0. The Macauley bracket «r w > has either a value of 0 whenever 
|S-J3 1 < y (defining the elastic or reversible domain) or a value of cr v whenever [S' -B | > Y (defining the ine- 
lastic or irreversible domain), with |S-B| = Y establishing the yield surface. Viscoplasticity (a rate dependent 
theory) admits states within the yield surface, on the yield surface, and outside the yield surface. Plasticity (a 
rate independent theory), on the other hand, only admits states within and on the yield surface; it docs not admit 
states outside the yield surface. The expression for the Zener parameter given above is very general, and 
includes three important special cases: i) viscoplastic theories without a yield surface (i.e., Y= 0) as used by 
MILLER [11], WALKER [20], KREMi’i. el al. [26], and others; u) viscoplastic theories with no evolving drag 
strength (i.e., D= constant ) as used by CHABOCliE [13], MARQUIS [18], and others; and tit) viscoplastic theories 
where the yield and drag strengths arc proportional (i.e., Y<*D) as used by PERZYNA [7], ROBINSON [15], and 
NOUAILHAS [31]. 

Phenomenologically, each internal variable /A fc {B tJ Jd ,Y } is taken to competitively evolve through a pro- 
cess of the type (cf. chaboche [13,39], walker [20], and miller et al. [27,38]) 

A - hardening - dynamic recovery - static recovery (29) 
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where the hardening term accounts for strengthening mechanisms, and the recovery terms account for softening 
mechanisms. The hardening and dynamic recovery terms both evolve with inelasticity; they are strain induced. 
The static recovery term evolves thermally; it is time induced. What functional form each particular term takes 
on varies with the type of internal variable, as shown below. 

At this time we list the basic structure for a thermodynamically admissible theory of viscoplasticity. A 
detailed discussion of the evolution equations follows in the next section. Our general theory of viscoplasticity 
for initially isotropic materials is characterized by the constitutive equations 


Sij = 2 Ufa; - e*) , CTii = 3 K (e„ - a A7" 8„ ) 

P.5 

5=i 5 5=i 

D - D 0 = £ D t = £ 6 ; 

5=i 5=i 

Y-To = i,Yj=I l hjy 1 

Y=1 T=1 

by the evolution equations 


ifj = 0 Z n,, - |e p l rty 



f Dr 


5^ = c 

l'-f 


yy = 

[-£] 

1^1 -Or/ 

with directions (unit or projection) defined as 



3 

n>1 - 2 |5-B | ’ 

b -5_ 3 B ‘> „5 = 

" 2 |B*| * " 

d,} - (\-p0bfj + P\ ^ n, 


2 


and therefore, by the intrinsic dissipation inequality 


(30) 


(31) 


(32) 
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v Dr rr jr Yy r y 


Here n = Et 2(l+v) is the elastic shear modulus; K = £73(l-2v) is the elastic bulk modulus; 
c,j = tij - (e u /3)5 i; is the deviatoric strain; (3,,, 5 and y are the thermodynamic fluxes conjugate to the ther- 
modynamic forces Bij, D-D 0 , and Y-Y 0 ; and H , h , and h' are their associated hardening coefficients. The 
material constants of the theory are: D 0 > 0 is the initial value of drag strength, Y o >0 is the initial value of 
yield strength, p^e [0,1] is a nonproportionality parameter associated with the constituent of back stress, 
and ce(0,l] is a parameter that can be arbitrarily set to ensure satisfaction of the intrinsic dissipation inequal- 
ity, equation (33). In the virgin or annealled state, - 0, D = D 0 , and T = To! therefore, B,j = 0, Dr = 0, 
and Y y = 0 in this state. The material functions of the theory are: 0 > 0 is the thermal difusivity; 

Z = |e p |/0 > 0 is the Zener parameter; L > 0, / > 0, and /' > 0 are the limiting stales for the internal vari- 
ables, as associated with dynamic recovery; and R > 0, r > 0, and r' > 0 are the static recovery functions. 
The limit and static-recovery functions are, in general, functions of state as established by the set 
{.S, ; J],j ,D ,Y ,T], while the functional dependences of thermal diffusivity and Zener’s parameter are specified in 
equations (27 and 28). These constitutive, evolution, and dissipation equations are a special case of the thermo- 
dynamic relationships given in equations (16-18). 

The kind and the number of independent state variables chosen, and the specific forms chosen for their 
associated material functions, vary from model to model. A particular model may be simple and employ only a 
couple facets of the theory, or it may be complex and employ many facets of the theory. This discretion is left 
to the individual modeler who best understands his own needs. 


In addition to equations (30-33) which characterize the deformation properties of the material, there arc 
the constitutive equation for entropy, the evolution equation for heat flux, and the thermal dissipation inequality 
given in equations (16-18) which characterize the thermal properties of the material. 


5. DISCUSSION OF THE THEORY 

An arbitrary number ( i.e u) of constituent back stresses = 2//^/3 can be introduced into a visco- 
plastic model. It has been demonstrated by CHABOCHE et al. [16,22,35,39], WALKER [20], MILLER et al. 
[27,29,38], and MCDOWELL & moosbrugger [30] that two back stresses (sometimes referred to as short-range 
and long-range back stresses) are typically necessary for proper material characterization. They arc thermo- 
dynamically admissible within this theoretical framework provided that i) the back stress is given by the entire 

U 

sum of its constituents = X s .) • and that ii) each constituent evolves in a manner compatible with equation 

5=i 
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(31). The back stress is therefore a sum of shape or weight functions, each being its own, generalized, state- 
variable coordinate. 

The implementation of the phenomenology of competing hardening and recovery mechanisms into an evo- 
lution equation for the back stress is complicated by the fact that the back stress is a tensorial quantity, and 
therefore has directional characteristics. Here the back stress is considered to harden linearly in the same direc- 
tion that the inelastic strain evolves in, as advocated by PRAGER [51]. Static recovery is taken to occur in a 
direction that opposes itself, in accordance with the Bailey-Orowan model of ponter & leckie [12]. And for 
the dynamic recovery of back stress, we adopt the concept of BURLET & cailletaud [57], as it presents a natural 
generalization to the popular model of Armstrong & Frederick [8]. Figure 1 presents a schematic for the 
directional characteristics of this phenomenology. It is simplified here for the case of one constituent of back 
stress for illustrative purposes. The Armstrong-Frederick model considers dynamic recovery to occur in a direc- 
tion that opposes itself (like static recovery), resulting in a formulation that is equivalent to a two-surface theory 
of plasticity with Mroz hardening, as shown in reference 18. The Burlet-Cailletaud model is equivalent to the 
Armstrong-Frederick model for proportional loading histories, but it has the capability for predicting a different 
response under nonproportional conditions, depending upon the value of the nonproportionality parameter Pt, 
found in equation (32). The one extreme of pz > = 0 gives dynamic recovery as advocated by Armstrong and 
Frederick. The other extreme of p^= 1 causes dynamic recovery for the constituent of back stress to occur 
in a direction —nfj. This direction does not oppose the direction of back stress; rather, it opposes that of inelas- 
tic strain rate (for the case of one constituent in back stress, as in the Burlet-Cailletaud model), which general- 
izes Prager’s, linear, kinematic, hardening rule into a nonlinear, kinematic, hardening rule. These two extremes 
can predict vastly different nonproportional responses. For example, tension-torsion ratchetting behavior need 
not shakedown for the case of p ^ = 0, while shakedown occurs almost immediately for the case of = 1 [57]. 
Reality is likely to be somewhere between these two extremes, but probably closer to Armstrong-Frederick (i.e. 
Pf. = 0) than nonlinear Prager (i.e. p ^ = 1) [32]. The coefficient n^B^I\B^\ to the unit normal n§ appearing in 
equation (32) is necessary for two reasons. First, it is required to obtain an equivalent predicted response under 
proportional loading conditions for all values of p ^ And second, it guarentees a positive-valued contribution to 
the intrinsic dissipation inequality (i.e., (n^B^j) 2 since n^B^j need not be - and often is not - positive valued). 

The basic form of the evolutionary equation for back stress, as proposed by Armstrong & Frederick [8], 
has been used by chaboche et al [13,16,22], WALKER [20], and many others in their viscoplastic models. In the 
absence of static recovery, this relationship introduces an evanescent strain memory effect (evanescent along the 
inelastic strain path) caused by a strain-induced dynamic recovery phenomenon which competes against the 
mechanism of linear strain hardening. This results in the back stress B^ evolving asymptotically to a limit 
state L associated with kinematic saturation. Considering just one constituent for illustration, one can combine 
the constitutive and evolution equations (30 and 31) for back stress into a form of the type proposed by 
Armstong and Frederick; it is 

B tJ = | \[efj - <(> \Ef\ by] (34) 
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under isothermal conditions. (Nonisothermal conditions would bring about an influence due to the temperature 
dependence of the hardening coefficient H - cf. CHABOCHE [13,39] and WALKER [20].) Convergence is 
exponential whenever X and <J> are constants [8], but in our general theory 


X = // 


1 -p 




(HO 




and <(> = 


R 

+ Z 


1 - p 




(35) 


Hence, the model of BURLET & cailletaud [57] effectively changes the concept of linear strain hardening into 
one of nonlinear strain hardening for the evolution of back stress. A graphic illustration of this evanescent con- 
cept is given in figure 2. Here the limiting state of back stress defines a hypersurfacc <}> = 1 in state space 
where | B | = B sa , . Whenever this bound is reached, a perfectly plastic or steady-state creep response is attained. 
In the absence of static recovery, B sal = L where L is the rate-independent (or plastic) limiting state of back 
stress; that is, L is the upper bound of B sa , . The value of this limit state is known to vary as the strength of 
the material varies [18,20], i.e., L(D,Y,T). The presence of static recovery reduces the limiting state of back 
stress making it rate-dependent, such that 0 < B sat < L. In the absence of dynamic recovery, the evolution equa- 
tion (31) is in accord with the BATLEY-OROWAN [58,59] hypothesis, that creep deformation evolves through □ 
competitive process between the mechanisms of strain hardening and time-induced static (or thermal) recovery. 
At kinematic saturation, the inelastic strain rate kjj becomes coaxial with both the deviatoric stress S tJ and the 
back stress B;j (in agreement with experiment [60,61]), and the nested set of flow surfaces 
( |5 —B | = constant } becomes stationary until unloading occurs. Otherwise, this set of flow surfaces can translate 
freely within the limiting hypersurface, as governed by the flow and evolution equations. 

Viscoplastic models use either inelastic strain (i.e. |c p I) or inelastic work (i.e. S tJ or (S,j ~B i; )?,, ) to 
drive the evolution of the isotropic variables D and Y. Note that (5, y -5 iy )e^ is the proper expression to use 
for inelastic work when a back stress is present, since in this case S,y£,/ is not always positive valued. A set 
of critical experiments has been performed by krempl [62] to determine which of these forms for the evolution 
of internal state is most representative of material response. He determined that inelastic strain provides the best 
characterization of strain hardening behavior. 

A hardening/recovery format for the evolution of drag strength D has been used by BODNER et al. [9,171, 
miller el al. [ 1 1,21,27,29], WALKER [20], and others in their viscoplastic models. Bodncr uses inelastic work to 
drive its evolution, while the others, for the most part, use inelastic strain for this purpose. Since a parameter 
ce(0,l] can be arbitrarily introduced into an evolution equation for drag strength (e.g. see equation (31)), cither 
hardening mechanism (i.e. inelastic strain or inelastic work) can be made to satisfy the intrinsic dissipation ine- 
quality of equation (11). Actually, the parameter c is transparent to the viscoplastic modeler, since he has the 
flexibility to arbitrarily adjust the hardening coefficient h found in the constitutive equation (30). There is 
only one driver for the evolution of drag strength that satisfies intrinsic dissipation independent of a parameter 
r; it is (Sij -Bij )ifjl(D -D 0 ) as introduced by FREED & WALKER [32], Multiple drag strengths are admissible 
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within the constructs of this theory, i.e. D = Dq + They have been advocated in the models of MILLER 

5=i 

et al. [21,27,29,38] to complement the short- and long-range back stresses, and to account for solute strengthen- 
ing effects. Except for the fact that the evolutionary equation for drag strength given in equations (30 and 31) is 
a scalar equation instead of a tensorial one, it has the same mathematical features as the evolutionary equation 
for back stress. 


The evolution equation for yield strength Y given in equations (30 and 31) is used in the models of cha 
BOCHE [13,16,22,39]. Our theory also permits the yield strength to be given as a sum of constituents 

w 

Y = Y 0 + however, there seems to be no physical significance in permitting w to exceed 1. The evolu- 

lion of yield strength has the same mathematical features as the evolution of drag strength, except that a parame- 
ter like c is not required to satisfy intrinsic dissipation. 

There are two types of viscoplastic models typically found in the literature; those with a yield surface, and 
those without one. Each is a special case of our theory. Each has a different energy-storage/heat-dissipalion 
diagram, as illustrated in figure 3. If one considers a material whose internal state is fully recovered or 
annealled, and if one deforms this material in monotonic tension in the absence of thermal recovery, then for 
those models with no evolving drag strength, i.e,, D = Dq, the intrinsic dissipation inequality (equation (33)) 
becomes 


v v B 2 < 7-Yo. f 

& + t 0 + -£~h p 


if >0 


while for those models with no yield surface, i.e., Y = Y 0 = 0, it becomes 


. B 2 (.d-d o ) 2 

c v - c (D -D 0 ) + — + C 


e p > 0 


(36) 


(37) 


Here only one constituent for each internal variable is considered for simplicity. These are the equations dep- 
icted in figures 3a and 3b, respectively. There the isotropic and kinematic variables are shown to saturate at the 
same level of accumulated inelastic strain, which is not an accurate representation of our observations, but it is a 
convenient one for illustrative purposes. These models predict energy-storage/heat-dissipation ratios that differ 
between themselves, but probably not to a large extent. As a consequence of effecting the intrinsic dissipation, 
the parameter c in the drag strength model influences the amount of inelastic work done on a body that is 
stored internally versus that which is dissipated away as heat. At saturation, both models predict that all the ine- 
lastic work done Si p is dissipated away as heat, in agreement with experiment [63]. 
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6. SUMMARY 


Constitutive, evolution, and dissipation equations governing the response of a general class of materials are 
derived from physical and thermodynamical considerations. These equations apply to materials that are initially 
isotropic, and whose strain is additively decomposed into elastic (reversible) and inelastic (irreversible) com- 
ponents. A general structure for a theory of viscoplasticity is then proposed. A second invariant formulation 
with three fundamental internal state variables is considered. One internal variable (the tensor-valued back 
stress) accounts for kinematic or flow-induced anisotropic effects. The other two variables (the scalar-valued 
drag and yield strengths) account for isotropic effects. Each can be composed as a sum of independently evolv- 
ing constituents. All the internal variables are assumed to evolve phenomenologically through a competitive pro- 
cess of strain hardening, strain-induced dynamic recovery, and time-induced static recovery. The thermodynamic 
constraint of intrinsic dissipation requires attention only for the hardening term of drag strength. By introducing 
an arbitrary parameter c into the evolution equation for drag strength, the constraint of intrinsic dissipation can 
be satisfied. 
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Figure 1 . - Schematic of directions for linear hardening LH , 
dynamic recovery DR, and thermal recovery TR of back 
stress. 
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INELASTIC STRAIN 

(b) Back stress - drag strength model (Y = O). 
Figure 3. - Energy storage and heat dissipation due to 
inelastic deformation. 
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